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Sidles and Sigg have shown that advanced LIGO interferometers will encounter a serious tilt 
instability, in which symmetric tilts of the mirrors of an arm cavity cause the cavity's light beam 
to slide sideways, so its radiation pressure exerts a torque that increases the tilt. Sidles and Sigg 
showed that the strength T of this torque is 26.2 times greater for advanced LIGO's baseline cavities 
— nearly flat spherical mirrors which support Gaussian beams ("FG" cavities), than for nearly 
concentric spherical mirrors which support Gaussian beams with the same diffraction losses as the 
baseline case — "CG" cavities: y^G/y 00 — 26.2. This has motivated a proposal to change the 
baseline design to nearly concentric, spherical mirrors. In order to reduce thermoelastic noise in 
advanced LIGO, O'Shaughnessy and Thorne have proposed replacing the spherical mirrors and 
their Gaussian beams by "Mexican-Hat" (MH) shaped mirrors which support flat-topped, "mesa" 
shaped beams. In this paper we compute the tilt-instability torque for advanced- LIGO cavities with 
nearly flat MH mirrors and mesa beams ( "FM" cavities) and nearly concentric MH mirrors and mesa 
beams ("CM" cavities), with the same diffraction losses as in the baseline FG case. We find that the 
relative sizes of the restoring torques are T CM /T CG = 0.91, T FM /T CG = 96, T FM /T FG = 3.67. 
Thus, the nearly concentric MH mirrors have a weaker tilt instability than any other configuration. 
Their thermoelastic noise is the same as for nearly flat MH mirrors, and is much lower than for 
spherical mirrors. 



I. INTRODUCTION 

Thermoelastic noise is the dominant noise source in 
advanced-LIGO interferometers with sapphire mirrors, 
at and somewhat below the frequency of optimal sen- 
sitivity 0, 0. O'Shaughnessy and Thorne H have 
proposed lowering this thermoelastic noise by flattening 
the cross-sectional profile of the arm cavities' light beams 
— i.e., by replacing the standard Gaussian-shaped beams 
by "mesa" -shaped beams (thick curves in Fig. [5] below). 
This can be achieved by replacing LIGO's nearly flat, 
spherically shaped mirrors by mirrors that have a nearly 
flat " Mexican- hat" (MH) shape. O'Shaughnessy, Strigin, 
and Vyatchanin have shown that the resulting re- 

duction of thermoelastic noise power is about a factor 3, 
and with d'Ambrosio and Thorne, they have carried out 
extensive studies which suggests that MH mirrors have 
no serious disadvantages compared to spherical mirrors 

Sinn. 

Sidles and Sigg0, @ have recently rediscovered a tilt 
instability in Fabry-Perot (FP) cavities, first pointed out 
by Braginsky and Manukin [|J , and they have shown that 
this instability is a serious issue for advanced LIGO's 
arm cavities, because of their high circulating light power 
(about 800 kW) and resulting high light pressure. In this 
instability, random forces cause the cavity's mirrors to 
tilt in a symmetric way 1 (Fig. and this tilt causes 



the light beam to slide sideways in the cavity by the dis- 
tance 8x sym shown in the figure, so its light pressure ex- 
erts a torque T on the mirrors that tries to increase their 
tilt. [Sidles and Sigg also showed that, when the mir- 
rors are tilted in an antisymmetric way as in Fig.^i, the 
resulting torque is stabilizing rather than destabilizing.] 
Sidles and Sigg analyzed the tilt instability, using geo- 
metric arguments, for cavities with nearly flat, spherical 
mirrors and their Gaussian light beams ("FG" cavities), 
and also for nearly concentric, spherical mirrors and their 
Gaussian beams ("CG" cavities). [The mirrors must be 
nearly flat or nearly concentric in order to make the light 
beams significantly larger than the Fresnel diffraction 
size, b = \J AL/2ir with A the light's wavelength and L 
the cavity length; large beams are required to keep the 
thermoelastic noise small.] Sidles and Sigg found that 
the instability is much more severe for the baseline FG 
cavities than for CG cavities with the same beam radii 
at the mirrors and thence the same diffraction losses. On 
this basis, the baseline design for advanced LIGO 01 has 
been changed from FG cavities with nearly flat mirrors 
to CG cavities with nearly concentric mirrors. 

Motivated by this Sidles-Sigg work, Thorne has pro- 
posed a mathematical way to design nearly concentric 
MH mirrors that support mesa beams with precisely the 
same mesa-shaped light-power distributions on the mir- 
rors as for the original nearly flat MH mirrors. Thome's 



1 Sidles and Sigg 0, use the opposite convention from us for 



"symmetric and "antisymmetric" tilt. 
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FIG. 1: FP resonator with (a) perfectly positioned spherical 
mirrors, (b) symmetrically tilted spherical mirrors and (c) 
antisymmetrically tilted spherical mirrors 



mathematical construction is presented, along with some 
generalizations of it, in a companion paper by Bon- 
darescu [TT| . 

In the present paper, we analyze the tilt instability 
for advanced-LIGO arm cavities with (i) nearly flat MH 
mirrors and their mesa beams ("FM" cavities), and (ii) 
with Thome's new nearly concentric MH mirrors and 
their mesa beams ("CM" cavities). We employ first- 
order perturbation theory in our analysis, by contrast 
with the Sidles-Sigg geometric techniques. We compare 
the strength of the tilt's destabilizing torque T for FG, 
FM, CG, and CM cavities that have beam sizes chosen so 
they all have the same diffraction losses, about 20 ppm; 
and we explore two choices for the radius of the mirror 
coating on the substrates: the baseline radius (14.7 cm), 
and a larger coated radius (16 cm) used in the analysis 
of d'Ambrosio et. al. 0>0| [their fiducial configuration] . 

In our numerical solutions to the eigenequation for the 
light's eigenmodes inside FM and CM cavities, we discov- 
ered remarkable duality relations between cavities with 
axisymmctric mirrors that deviate by an amount H (r) 
from flatness, and cavities with mirrors that deviate by 
—H(r) from concentric spheres. We verified these nu- 
merically discovered duality relations for several different 
forms of H(r), in addition to those of MH mirrors. This 
motivated Chen and Savov, and independently Agresti 



and d'Ambrosio ^| to devise analytic proofs of our du- 
ality relations. The duality relations provide a unique 
one-to-one mapping between the eigenstates and eigen- 
values of the dual cavities — a mapping that may be 
useful not only for advanced LIGO but in a variety of 
other applications of Fabry Perot cavities. 

This paper is organized as follows. In Sec. II, we use a 
first-order modal analysis of a Fabry-Perot cavity to de- 
rive a general formula for the torque exerted on the mir- 
rors when the cavity is perturbed, in terms of as-yet un- 
known mode coupling coefficients o-k and mode-overlap 
integrals Ik- In Sec. Ill, we use first-order perturbation 
theory of Gaussian-beam (FG and CG) cavities to de- 
rive analytical formulas for ctk and Ik, and thence for 
the tilt-induced torque T in the FG and CG cases, and 
we show that our formula for the torque is equivalent to 
that of Sidles and Sigg 0,0 In Sec. IV we use first-order 
perturbation theory to derive formulas for the coupling 
coefficients a&, and thence for the torque T, in terms of 
a cavity's eigenvalues and mode-overlap integrals Ik- In 
Sec. V we present our numerical results for the modes 
and their eigenvalues for FM, CM, FG, and CG cavities, 
and we discuss the duality relations between the nearly 
flat and nearly concentric cases. Finally, in Sec. VI, we 
combine the numerical results of Sec. V with the formu- 
las of Sees. Ill and IV, to deduce the tilt-induced torque 
for our four cavity designs — using two sets of param- 
eters: those for cavities with advanced-LIGO baseline 
mirror radii, and those for d'Ambrosio et. al.'s slightly 
larger mirrors ("fiducial" configurations). We present a 
brief conclusion in Sec. VII. For the readers interested in 
our numerical implementation of the eigenvalue problem, 
we include an Appendix where we sketch details of our 
computational work. 

The results presented in this paper are based on previ- 
ous work on nearly flat configurations by S. Vyatchanin 
|13| (some errors in this paper are corrected here) com- 
bined with recent analyses of nearly concentric cavities 
by P. Savov. An analytical proof of the duality relation 
between nearly flat and nearly confocal resonators by P. 
Savov and Y. Chen, and independently by E. D'Ambrosio 
and J. Agresti, will be provided in a companion paper 

m 



II. MAIN FORMULAS 

The light inside LIGO arm cavities is well-described by 
the laws of diffraction optics in the paraxial approxima- 
tion. The eigenvalue problem in this approximation for a 
half trip through a cavity with two identical axisymmet- 
ric mirrors can be written as 



G(ri,r 2 )u(r 2 )d 2 r 2 = Aw(ri). 



(2.1) 



In the above equation u(r) is an eigenmode of the cavity 
and A is the corresponding eigenvalue. The eigenmode 
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represents the state of the light (the electric field) on the 
surface of a mirror. 

For advanced LIGO diffraction losses will be very small 
(about 10 ppm for each half trip), so it is an excellent ide- 
alization to ignore the losses and idealize the mirrors as 
infinite in radius. Then, |A| = 1, G is a unitary operator, 
and its eigenvectors form a complete set. Each eigen- 
mode u nm and the corresponding eigenvalue X nm are la- 
beled by two (quantum) numbers — radial (or principle) 
number n = 0, 1, . . . and angular (or azimuthal) number 
to = 0, 1, . . .. All modes with angular number to = 
are axisymmetric (no angular dependence), m = 1 are 
dipolar, m — 2 are quadrupole, etc.: 



oc e 



(2.2) 



The eigenmodes are normalized and orthogonal to each 
other according to the following definition: 



i (r) d T — ^mri2 ^mima 



(2.3) 



We will use this set of eigenvectors as a basis for ex- 
panding the eigenmodes of cavities with tilted mirrors. 
The radial coordinate r is dimensionless and measured 
in units of the Fresnel diffraction size 



b = y/L\/2lT. 



(2.4) 



When the mirrors of a FP cavity are tilted in a sym- 
metric way (as in Fig.^i), the cavity's fundamental mode 
Moo(^) is transformed into the fundamental mode uoo(r) 
of the perturbed cavity. The torque acting on the mirrors 
when the light is in this mode and has power P is 



T 



2Pb 



I Uqq (r) I r cos ipd 2 f. 



(2.5) 



The new fundamental mode can be expanded over the 
set of orthonormal modes {u nm (r}} of the unperturbed 
cavity 



u o(r) = u 00 (r) + 



(2.6) 



In this paper, we study effects only to first order in the 
perturbation. That is why the coefficient in front of 
uoo(f), in Eq. I|2.6[) . is unity. 

By substituting Eq. (|2.6[) into Eq. (|2.5|) and using the 
angular dependence of the eigenmodes Eq. I|2.2|l . we con- 
clude that only the dipolar eigenmodes (to = 1) con- 
tribute to the net torque and more specifically their part 
proportional to cos(y). Thus, for our purposes of calcu- 
lating the torque, we will assume u n \ oc cos(y). Since 
the only modes we use from now on are the fundamental 
mode uoo and all dipolar modes u n i, in order to simplify 
notation, we collapse the indices into one labeling index 

k = n + m. (2.7) 



Thus the fundamental mode becomes uo, the first dipolar 
mode becomes u\ (corresponding to the old notation moi) 
and so on. When necessary, we will use the conventional 
notation with two labeling indices. 

We will study the effects of tilt only to first order in 
the tilt angle 9, so for our purposes we use the following 
expansion of the perturbed eigenmode: 



fc=i 



uo{r) 



u (r) = u a (r) +^a k u k (f), (2 

Uk{r) 
\uo(r)] 2 rdr 



V2n 

u k (r) cos if 



[u k (r)] 2 rdr 



(2.9) 
(2.10) 
(2.11) 



fc=l,2,. 



(2.12) 



In the above equations, u k (r) are the dipolar modes on 
the surface of a mirror; u k (r) are their parts depending 
only on the radial coordinate r; all u k (r) arc dimension- 
less and normalized as shown above [cf. Eq. I|2.3|l with 
to = 0, l,k = n + m]; and a k are dimensionless coupling 
constants, proportional to the mirrors' tilt angle 9, which 
we will evaluate in Sec. Ill for Gaussian (FG and CG) 
beams and in Sec. IV for mesa (FM and CM) beams. In 
general, u k {r) are complex fields, but since the mirror 
surfaces coincide with the beam's wave front, up to an 
overall complex phase which we chose to be zero, they 
are real fields. 

Now we can calculate the torque that the cavity's light 
exerts on each mirror: 

T = / "[wo(r)] 5 V cos ip r dr dip — 

= 2^^Q!fe / uo(r)ufe(r) r cos ipr dr dip , 

C fc=i ^ 

where we have used Eq. I|2.8(l . By inserting Eqs. Ij2.9|l 
and H2.10f> . we obtain the following formulas for the tilt- 
induced torque to first order in a k (first order in 9): 



2V2Pb, ^ 
T = x y a k I k 

c 

k=l 



h. 



uo(r)u k (r) r 2 dr . 



(2.13) 
(2.14) 



These formulas are valid for any FP cavity and in par- 
ticular for FG, FM, CG, and CM cavities that interest 
us (of course, the modes u k (r) are different for different 
cavities) . 

In the sections below, we calculate the values for the 
coupling constants a k and the overlap integrals I k for 
our four types of cavities. Our analysis for conventional 
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spherical mirrors (FG and CG; Sec. Ill) is entirely an- 
alytical, whereas for any generic mirror shape, and MH 
mirrors in particular (FM and CM; Sec IV), numerical 
treatment is required. We will test our numerical meth- 
ods by applying them to FG and CG cavities and com- 
paring with the analytical results. 



Sx s 



III. GAUSSIAN-BEAM (FG AND CG) 
CAVITIES 



We consider a cavity with identical spherical mirrors. 
We are interested in a symmetric tilt of the two mirrors 
by a small angle 9 as shown on Fig. In this case, 
the axis of the new mode Uo(r) is displaced by a small 
distance 5x synl , but is still parallel to the old axis. The 
field distribution on each mirror will be unchanged, but 
shifted by Sx sym . 

Spherical cavities have been studied thoroughly (see 



spheric 

.14]); 



e.g I: their fundamental modes are the well-known 
Gauss-Laguerre modes (called in this paper FG and CG 
modes). We will use these modes derive analytical for- 
mulas for ak and Ik- The main axisymmetric and dipolar 
modes [uo(r) and Ui(r)] are given by (see e.g. (lif): 



uf(r) 



r = 



(1-02)1/4- 



(3.1) 
(3.2) 
(3.3) 



Here r is the dimensionless radial coordinate (mea- 
sured in units of b), r is the dimensionless radius of 
the beam at the mirrors' surface (also in units of 6), 
g = 1 — L/R is the so-called (^-parameter of the cav- 
ity, L is the distance between the mirrors, and R is the 

mirrors' radius of curvature (Fig.^). (The intensity on 

2 / 2 

the mirror is proportional to e~ r /r o.) 

For spherical mirrors the displacement of the optic axis 
(see Fig.Hb): 



RB _ L9 

sym ~ b ~ b(i- g y 



(3.4) 



Next, we write down the main mode Uq of the FP 
resonator with tilted mirrors and expand it to first order 



= ( r cos (p — So 



syni 



2 • 2 

r sin <p, 



(3.5) 
(3.6) 



l r 5x sym cosip \ 



Gr-* , ^sym COS (f q . . 



wo (r) 



2wr 

Sx sym uf(r) cos tp 
y/2r a/tt 



(3.7) 



As we can see, the only nonzero coupling constant is af 
L 9{1 +g) 1 /* 



V2b{l-gf/ A ' 



(3.8) 



From Eqs. $ZTfy . JOJ, (|3~2l . and @, we can easily 
calculate the only overlap integral we need for Gaussian 
beams: 



!? = / u o( r ) u i( r ) r2 dr 
Jo 

1 

" r " - (l-.g2)i/4- 



(3.9) 



Substituting into Eq. (|2.13|l along with Eq. (|3.8I) we de- 
rive a final expression for the torque: 



2pl e 

~ (1 - .9) ' 



(3.10) 



This result, derived by a modal analysis, is in complete 
agreement with the result of the Sidles-Sigg geometrical 
analysis in its long-cavity limit (Section 5 of In their 
notation, the torque for the unstable configuration is 



-fe_0 



2PL 9 



(1-9) 



(3.11) 



where = fc_ is the negative eigenvalue of a torsional stiff- 
ness matrix (Eq. (23) of Section 5 in jg). (Note that 
negative eigenvalues in the Sidles-Sigg analysis are as- 
sociated with unstable configurations — the subject of 
interest in this paper.) Our perturbation method gives 
the exact result (to first order in 9) for spherical mirrors, 
because the only contribution to the torque is from the 
lowest dipolar mode u\. This is a property only for spher- 
ical mirrors and their Gaussian beams. As we'll see in 
the following sections, for any generic mirror shapes, we 
have to calculate the contribution from all higher dipolar 
modes. 
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IV. MESA-BEAM (FM AND CM) CAVITIES: 
ANALYTICAL FORMULAS 

a. Perfectly positioned mirrors (Fig. QJiJ. For any 
cavity with axisymmetric mirrors, and in particular MH 
mirrors, the main axisymmetric mode uo(r) and all dipo- 
lar modes Ufc(r) satisfy the integral eigenequations 



G(fi, r 2 ) u a (r 2 ) d 2 f 2 = AoUo(n), (4.1) 
G(fi,f 2 )wfc(ri)rf 2 ri = A fe u fc (r 2 ), (4.2) 



where G, uo,Uk,r{, r 2 are all dimensionless and the eigen- 
value of the k th dipolar mode u k is 

In the paraxial approximation, the kernel of the oper- 
ator G is the following (up to a trivial factor of e lkL due 
to phase accumulation along the arm length L, which we 
omit, thereby fixing a common overall phase factor in all 
the A fc ): 



G(ri,r 2 ) = — exp 

Z7T 



h xa {r) = kH h2 (r), k = — . 



' I hi[ri) - hi{r 2 ) 

2tt 



(4.3) 



Here Hi(fi) and H 2 (r 2 ) are the physical deviations of the 
mirrors' surfaces from a plane surface, which we assume 
to be the same, Hi(r\) — H 2 (f 2 ) (identical mirrors). 

b. Symmetrically tilted mirrors (Fig. \Tfr). The tilt 
is equivalent to small deviations of each mirror's position 
from the unperturbed one: 

Shi = kbr\ cosy?! 9 (left mirror) (4.4a) 
5h 2 = kbr 2 cos(p 2 8 (right mirror). (4.4b) 

These tilts induce a coupling of all the dipolar modes 
ui,u 2l ... into the cavity's fundamental mode uq, as 
shown in Eq. (|2.8[) . though (as our numerical work will 
show) the coupling for the first dipolar mode is far greater 
than the others ol\Ii S> oi k l k for k = 2, 3, . . .. 

For simplicity, we will show the analysis only for the 
first dipolar mode u% (uqi in the conventional notation). 
The generalization for the higher dipolar modes is triv- 
ially obtained by by replacing the subscript 1 by the de- 
sired dipolar mode's subscript k 

The eigenvalue of the fundamental mode of the per- 
turbed cavity Ao will slightly differ from Ao: Ao = Ao + A. 
Thus, we have the following integral eigenequation for 
Uo(r) 

(A + A)[u (n) + axuiin)} = 

= J G(ri,r 2 )[l-i^ 1 (ri)-W/ l2 (f 2 )] 

x [uo(f 2 ) + aiui(r 2 )} d 2 r 2 . (4.5) 
This equation can be simplified by use of the eigenequa- 



tion of the original unperturbed system (|4.1() : 

Auo(r'i) + (A + A - X^axUxifx) = 

= -if G(r 1 ,r 2 )[^ 1 (r 1 ) + ^ 2 (r 2 )] 

x[u (r) +a 1 ui(f)]d 2 f 2 . (4.6) 

Multiplying Eq. (|4.6|l by ito(ri) and integrating over d n., 
one can find that the correction A to the eigenvalue Ao 
has second order of smallness, A ~ 9 2 , so below we set 
A = 0. 

Multiplying Eq. (|4.(i|l by ui(r*i) and integrating over 
d 2 r*i, one can find a\\ 

(A -Ai)ai = -i(A + Ai) 

x I u^r^uiir^Sh^r^d 2 ^, 



so 



Oil = 



ikbO (A + Ai 
^(Aq — Ai) Jo 



Uo(r) ui(r) r 2 dr 



^Ly^A^fAi) 
\/26(Ao-Ai) ' 

Similarly for the higher dipolar modes 



iLI k 6{\ + \ k ) 
V2b(X -X k ) 



(4.7) 



(4.8) 



In order to calculate the numerical value of a k , we 
must solve the eigenequations l|4.1|l and l|4.2|l numeri- 
cally for the eigenvalues Ao, Afc and the corresponding 
eigenfunctions uo(r), u k (r) (see Appendix B for details). 
The value of the integral I k can be calculated numerically 
from Eq. H2.14jl . 

Note that the formulas in this section are valid for 
any resonators with symmetric mirrors iii(ri) = H 2 {r 2 ) 
and very low diffraction losses, not just for mesa-beam 
resonators. 



V. NUMERICAL SOLUTIONS OF 
EIGENEQUATIONS 

We have solved the eigenequations (|4.1|l and (|4.2I) nu- 
merically using the scheme described in Appendix B for 
our four cavity configurations: FG, CG, FM, and CM. 
Recall that our nearly flat and nearly concentric cavities 
were chosen such that the intensity u$(r) 2 and therefore 
uo(r) at the mirrors' surfaces are identical (FG and CG 
are the same and FM and CM are the same). We have 
found numerically for FM and CM (mesa beams) and 
for FG and CG (Gaussian beams, Sec. Ill) that Ufc(r) is 
also the same for the nearly concentric and nearly flat 
cases. The eigenfunctions uq and u± are shown in Fig. [21 
The eigenvalues, by contrast, are different for nearly flat 
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FIG. 2: Fundamental modes uo(r) (thick curves) and first 
dipolar modes iti(r) (thin curves) at mirrors' surfaces for (a) 
FG and CG cavities, and (b) FM and CM cavities. The modes 
are dimensionless and normalized according to Eqs. 12. lit . 
and l2.12H . We have used the fiducial cavity parameters of 
d'Ambrosio et. al.: Eqs. (2) of Sec. IVA of and Sec. IIIA 
of H 



and nearly concentric cavities, so we have four sets of 
eigenvalues (FG, CG, FM, CM), depicted in Fig. El 

In our numerical solutions to the eigenequations (|4.1|) 
and 1)4. 2|) . we found an interesting duality relation be- 
tween nearly flat and nearly concentric configurations. 
This duality relation is satisfied for any generic mir- 
ror shape that satisfies the paraxial approximation. To 
within numerical error of less than 0.05 per cent, we 
found that a nearly concentric cavity, which has the same 
intensity profile as a nearly flat configuration, also has 
the same mirror-shape correction as the nearly flat cav- 
ity, but with opposite sign: 



8h c (r) = -5h F (r). 



(5-1) 



Here Sh c (r) is the deviation from concentric spherical 
shape, and Sh F (r) is the deviation from flat shape. We 
also found, numerically, a unique mapping between the 
eigenvalues of these dual configurations: 




FIG. 3: Eigenvalue spectrum in the complex plain. Note 
that all eigenvalues satisfy the duality relation, Eq. ]5.2H (n — 
0, m = for Ao, and n — 0, m = 1 for Ai); see also Hal . 



for any pair of integers n,m — 0, 1, 2, . . .. In addition, 
all higher modes have the same intensity profiles at the 
mirrors' surfaces as their counterparts 

\u C nm \ 2 = \< m \ 2 (5.3) 

for any integer n, m = 0, 1, . . .. 

Remarkably, our numerical calculations showed that 
these relations hold not just for mesa-beam cavities, but 
for all stable cavities that we explored (all mirror shapes 
Sh C F , including cavities in which the deviations 5h G ' F 
from concentric spherical and flat shapes are large — so 
large as the paraxial approximation is valid) . 

This has led us to conjecture a duality relation between 
symmetric cavities with axisymmetric mirrors: for any 
two such cavities, A and B, with 



h A {r) + h B (r) = — 

Li 



(5.4) 



there exists a one-to-one correspondence between their 
eigenstates: they all have the same intensity profiles at 
the mirrors, while 



Kim — ( - l) m+1 (^nm) 



(5.5) 



(5.2) 



Chen and Savov, and independently Agresti and 
d'Ambrosio [12| have verified this conjecture analytically. 



VI. STRENGTH OF THE TILT INSTABILITY 
FOR FG, CG, FM, AND CM CAVITIES 

We now have all the tools we need to compute the tilt- 
induced torque T on the cavity's mirrors, for FG, CG, 
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TABLE I: Comparison Between Analytical and Numerical Re- TABLE II: Numerical Results for FM and CM cavities; 
suits for FG and CG Cavities; Qi is measured in units of a k is measured in units of (#/10 -8 ) and T& is in units of 
(<9/l(T 8 ) and T is in units of (Pb/c)(6/W- a ) (Pb/c)(6/W- s ) 





Analytical 


Numerical 


k 


T FM 

i k 


„ FM 

a k 


rpFM 
1 k 


rCM 


„CM 

a k 


rpCM 




FG 


CG 


FG 


CG 


1 


2.6464 


0.04525 


0.33867 


2.6464 


0.00018 


0.00137 


h 


1.8075 


1.8075 


1.8073 


1.8073 


2 


0.1136 


0.00009 


0.00003 


0.1136 


0.00016 


0.00005 


Q.1 


0.012526 


0.00030802 


0.012525 


0.00030799 


3 


-0.015 


-0.00000 


0.00000 


-0.015 


-0.0002 


0.00001 


T 


0.064038 


0.0015747 


0.064023 


0.0015743 


Total 






0.33870 






0.00143 



FM, and CM configurations. We shall evaluate T for two 
sets of cavity parameters: the fiducial parameters used 
by d'Ambrosio et. al. [J 0] and the advanced-LIGO 
baseline parameters (Table 1 in 

The set of parameters for the fiducial cavity (see 
Sec. IVA(2) of 3] and Sec. Ill A of 3) is: 

L — 4 km — the length of the cavity. 

A = 1064 nm — the wavelength of the laser beam. 

k = 2n/X — the wave number associated with A. 

b = ^LX/2n = 2.603 cm 
length scale (Fresnel length). 



the natural diffraction 



= 16 cm 



the radius of the mirrors' coated 



surfaces. 

g FG — 0.952 — the g-factor for the fiducial FG res- 
onator (corresponding mirror radius of curvature R — 
83.33 km). 

g CG = -0.952 — the ^-factor for the fiducial CG res- 
onator (corresponding mirror radius of curvature R = 
2.05 km). 

r = 6/(1 - ,g 2 ) 1/4 = 4.7 cm — the radius of the FG 
and CG beams at the mirrors. 

D = 46 = 10.4 cm — the radius parameter of the 
FM and CM beams at the mirrors (see Sec. II A and 
Sec. IVA(2) of H). 

The above beam radii were chosen so as to make the 
diffraction losses be about 20 parts per million (ppm). 
More specifically, they are 23 ppm for the FG and CG 
beams and 19 ppm for the FM and CM beams. 2 

From Eqs. (|3.8|) . I|3.9|l . and l|2.13|l we can calculate the 
integral I G , the coupling constant a G , and the torque 
T G for the FG and CG cavities. Our results are shown 
in the second and third column of Table |U 

We have already established an agreement between our 
analytically derived results using the modal analysis de- 
scribed in Section III and the Sidles-Sigg results derived 
from geometric considerations @. We can also test the 
numerical first-order perturbation methods that we de- 
veloped for arbitrary mirror shapes by applying them to 
our FG and CG cavities. By substituting our numerical 



results for u FG = u UG , uf G = u GG , X FG , A GG , Af G , and 



,FG — -.CG n.FG — 
~ . ' 1 _ 

A GG into Eqs. ijSlJJl, 63, an d if2~T3j) . we calculate the 
results shown in the last two columns of Table HJ These 



numerical results all agree with our analytical results to 
within 0.05 per cent, thus validating our numerical meth- 
ods. 

As was found by Sidles and Sigg, the CG configuration 
is significantly less unstable than its nearly flat counter- 
part FG. The analytical analysis (first two columns in 
Table [I] predicts 



T FG 1 + g 



FG 



R 



FG 



T GG 1 - g 



FG 



R GG 



40.667, 



(6.1) 



which is in agreement with the numerical result 40.667 
(last two columns). 

From the modal analysis applied to FG and CG cav- 
ities [Eqs. (|4.7(l and l|2.13[) ]. we deduce that, aside from 
factors that are the same for FG and CG, 



T G oc i 



A G 



= cot an 



(6.2) 



Here, the equality holds because |A| = 1 for all modes 
(negligible diffraction losses) and G X is the phase sepa- 
ration between A G and A G , i.e. the argument of A G /A G 
(in Fig. |3 we show (f)™ 1 for the FM cavity). Thus, 
Eq. 16.2f> is governed by the phase separation of the eigen- 
values A G and A G . As Fig. |21 shows the two eigenval- 
ues for the FG configuration are very close to each other 
so cotan(0Q G /2) ^> 1, whereas the phase separation of 
the eigenvalues for the CG configuration is close to n so 
cotan(0 G 1 G /2) < 1. This explains why T FG > T GG . 

Similarly to the above Gaussian analysis, we use our 
numerical results to compute the torques T FM and T GM 
for FM and CM cavities respectively. In this case, we 
must include the contributions from higher order dipo- 
lar modes (ui,i*2, and u 3 ). From Eqs. (|2.13t . I2.14|l . 
and (|4.8|l . we have calculated the integrals 



FM 



tCM 
1 k J 



the coupling constants ajl , , and the torques 
T FM ,T CM for the FM and CM cavities. Our results 
are shown in Table |nj Note that the dominant contri- 
bution to the torque comes from the first dipolar mode, 
k = 1; the higher modes give contributions of only a few 
per cent, at most. 

For mesa-beam resonators, as in the case of Gaussian- 
beam resonators, the nearly flat configuration (FM) is 
far more unstable than its nearly concentric counterpart 
(CM) 



2 We have deduced these diffraction losses from our numerical so- 
lutions of the cavity's eigenequation. 



(6.3) 
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TABLE III: Comparison between different configurations of 
a fiducial optical cavity. The torques due to light pressure 
(when tilt angle 9 and circulating power P are the same) are 
normalized such that T CG — 1. 





Nearly Flat Cavity 


Nearly Concentric Cavity 


G-Beam 


T" u = 40.7 




M-Beam 


T FM = 


T CM = 0.91 



In this case, the discrepancy is even bigger than in the 
Gaussian case since the eigenvalues for the FM configu- 
ration are closer to each other on the unit circle (Fig. 
than for the FG configuration and the phase separation 
of the eigenvalues for the CM configuration is even closer 
to 7r than the phase separation of the eigenvalues for the 
CG configuration (Fig. |3J. 

In Table IIIII we compare all four configurations FG, 
CG, FM, and CM, normalized by T CG . For nearly flat 
resonators, going from a Gaussian-beam to a mesa-beam 
configuration increases the strength of the instability by 
about a factor 5. There are two effects contributing to 
this increase as we can see from the following relation (in 
which we focus on the dominant, k = 1 contribution): 



af I M 
af I M 



A I 



A G + Af 
cotan (0oi/2 



cotan Uqi/2) 



A? 



I G 



Af 



M 



I G 



(6.4) 



In the case of the nearly flat configurations both phase 
differences are small and since (j)™ < <j) FG (see Fig.|3J), 



cotan (^ M /2) 
72) 



cotan I 



> 1. 



This effect is amplified by the second ratio because of 
the higher overlap between the two eigenstates in the 
case of mesa beams than for Gaussian beams. This is 
manifested in the higher value of I F =2.65 compared 
to I FG = 1.87 (compare the overlaps between each pair 
of modes in Fig. an d Fig. \^p). 

For nearly concentric resonators, going from Gaussian- 
beam resonators to mesa-beam resonators weakens the 
net instability: T CM /T CG = 0.91. In this case, the dif- 
ference in the overlaps of the eigenstates is unchanged, 
but the phase differences are close to, but less than ir. 
Since cj) ^ 1 > cj) GG (again look at the separation of each 
set of eigenvalues on the unit circle for the CG and CM 
configurations in Fig. [3, 



for nearly flat and nearly concentric cavities is straight- 
forward using Eg. I| 6 |.4 |l and the duality relation (see 
Eq. JS3 and Ref. [l2]). 

In our formulation of the perturbation theory, we ac- 
count for effects scaled to first order in the tilt angle 
6. We assume small mode mixing <C 1 in order for 
the perturbation method to work. From our numerical 
results (Table 111}) . we see that at <C I requires the an- 
gular orientation of the cavity mirrors be controlled to 
9 < Iff- 8 . 3 

The contributions Tj~ of the higher order dipolar modes 
k = 2, 3, . . . to the torque can be understood by studying 
the analog of Eq. (|6.4|) . From the relative locations of the 
eigenvalues along the unit circle and the overlapping of 
the eigenmodes, it is easy to show that 7Vs are monoton- 

ically decreasing, T\ > Ti > T3 Thus, we accept the 

contribution from the highest dipolar mode 1*3 in our cal- 
culation, including the numerical error, as the maximum 
error of the method due to neglecting the higher order 
dipolar modes. In this way, we conclude that the error in 
our total torque in the case of the CM cavity is less than 
I per cent. In the case of the FM cavity the error of the 
method is practically of order of the numerical error, so 
it is less than 0.1 per cent. 

For another comparison, we perform the same calcula- 
tions for the baseline design of advanced LIGO (Table I 
in |10l|). The baseline parameters were chosen such that 
the beam radius at the mirrors 4 in the case of spheri- 
cal mirrors is 4.24 cm, corresponding to diffraction losses 
of 10 ppm. The MH-mirror configurations are designed 
to have about the same diffraction losses. The resulting 
baseline parameters are: 

r max — (15.7 — 0.8) cm = 14.9 cm — the radius of the 
coated mirrors' surfaces. 

g FG — 0.9265 — the g-factor for FG resonator (corre- 
sponding mirror radius of curvature R — 54.44 km). 

g CG — —0.9265 — the g-factor for CG resonator (cor- 
responding mirror radius of curvature R = 2.076 km). 

r = b/(l - ,g 2 ) 1/4 = 4.24 cm — the radius of the 
Gaussian beam at the mirrors. 

D = 3.36 = 8.58 cm — the radius parameter of the 
mesa beam at the mirrors. 



Tabic IIVI contains the final results for these baseline 
parameters (including the sum of the contributions to 
the torques from the first three dipolar modes). Again, 
the least unstable configuration, and thus the easiest to 
control against tilt, is the nearly concentric mesa-beam 
(CM) resonator. 



cotan 



(<^ M /2) 



cotan I 



/2 



< 1. 



The two effects counteract each other and for this choice 
of parameters the net result is in favor of the CM- 
Beam resonator. The comparison between the torques 



3 Currently, the control system of the initial LIGO interferome- 
ters operates with accuracy 8 ~ 10~ 7 ; an a ccur acy ~ 10 — 8 is 
planned for advanced LIGO interferometers I Hi . 

4 Note tha t our definition for the beam radius at the mirrors differs 
from 111! by factor of y/2. 
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TABLE IV: Comparison between different configurations of 
a cavity with parameters of the current baseline design for 
advanced LIGO. The torques due to light pressure (when tilt 
angle 8 and circulating power P are the same) are normalized 
such that T CG = 1. 





Nearly Flat Cavity 


Nearly Concentric Cavity 


G-Beam 


T*' u = 26.2 




M-Beam 


T FM = gg 


T CM = q 



VII. CONCLUSIONS 

As Table ITvl shows . by switching from the conventional 
Gaussian-beam cavities for advanced LIGO to concentric 
mesa-beam cavities, the instability to symmetric tilt will 
be reduced (dramatically compared to a FG cavity and 
moderately compared to a CG cavity). Furthermore the 
sensitivity of the interferometer will improve significantly 
due to the reduced thermoelastic noise (See Table I in Q 
and also [10). 
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APPENDIX: NUMERICAL SOLUTIONS OF 
CAVITY EIGENEQUATIONS 

In order to generate the set of basis solutions needed 
to construct perturbation theory for a cavity with arbi- 
trary mirror shapes, we must numerically solve an in- 
tegral eigenequation. We have done so using the fol- 



lowing method, based on earlier work by O'Shaughnessy 
(Sec. VB of @). 

Since the mirrors are axisymmetric [h(r) = h{r)\, we 
can decouple the angular and radial dependences in the 
eigenequations. In the numerical implementation of the 
eigensolver we used the following definition: 

u nm {r) = u nm (r)e- m ^, m = 0,1,2,.... (A.l) 
Note that, for m = 0,1, this definition of the funda- 
mental radial mode uo{r) and the dipolar radial mode 
u\(r) differ from the definitions in Eqs. (|2.9|l and l|2.10|l . 
However, after solving the eigenequations, all modes are 
renormalized by numerically computing the integrals in 
Eqs. (|2.11|l and l|2.12|l . so at the end we have radial modes 
defined as in Eqs. I|2.9|l and l|2.10|l . The resulting Uk are 
the radial modes we need for computing Ik in Eq. I|2.14|l . 
By plugging Eq. HA.1|) into Eqs. I|4.1|l and (|4.2|l and in- 
tegrating over the azimuthal angle, we can reduce the 
eigenproblem to a one-dimensional integral equation 



A„ m u„ m (ri) = J Gm(r 1 ,r2)u nm (r2)r 2 dr2, (A. 2) 
G m (r u r 2 ) = (-i) m+1 J m (nr 2 ) 

x „p[i(£i±^-2*(r))" , 

where J m is the Bessel function of the first kind and order 
m. 

We discretize space along the mirrors' radial direction 
in a uniform grid 

rj =jr max /(N-l), j = 0,1,...,Y-1. (A.3) 

We define the matrix G( m )y = G m (xi, Xj), the eigen- 
vectors U(„jj = u n (xj) (m, n label the mode and i, j are 
indices to access the matrices' and vectors' components), 
and we approximate the integration by a simple quadra- 
ture rule. Then the integral eigenproblem reduces to a 
matrix eigenvalue problem: 



Anm u n = M m u ni with Af( m )y = — — j-G( m )jj. 

(A.4) 

This equation can be solved for X nm and u n by any stan- 
dard matrix eigensolution software package. 
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